The existence of a two-solar mass neutron star 
constrains the gravitational constant Gn at strong field 
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In General Relativity there is a maximum mass allowed for neutron stars that, if exceeded, entails 
their collapse into a black hole. Its precise value depends on details of the nuclear matter equation 
of state about which we are much more certain thanks to recent progress in low-energy effective 
theories. The discovery of a two-solar mass neutron star, near that maximum mass, when analyzed 
with modern equations of state, implies that Newton's gravitational constant in the star cannot 
exceed its value on Earth by more than 8% at 95% confidence level. This is a remarkable leap of ten 
orders of magnitude in the gravitational field intensity at which the constant has been constrained. 



The gravitational attraction force between two bodies 
of mass TTii and 7712 at distance r is, if not much precision 
be required, given by the renowned law of Newton 
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The Newtonian constant Gn, first measured by 
Cavendish, also features in the more precise field equa- 
tions of Einstein's General Relativity 
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that need to be used under intense gravitational fields 
with Einstein's curvature tensor i?^^ — -^RSi-lv and mat- 
ter source T^^ or otherwise if high precision is ex- 
pected. This constant has been carefully measured on 
numerous occasions [l|, |2[ and is currently taken to be 
6.6738(8)iV(m/fc5)2. 

The study of the orbital evolution of binary pulsars [3] 
has also allowed to establish the validity of General Rel- 
ativity (and indirectly, provided a measurement of Gn) 
under stronger gravity conditions. While on the Earth's 
surface g ~ 9.8 mj s^ , for the binary pulsar J0737-3039 
the relevant acceleration is g ~ 270 m/ s^ , and for PSR 
B1913-I-16, 330 m/s^. The assessment of Gm in those 
systems has a respective precision of 0.05% and 0.2% 
(from measurements of corrections to Kepler's law, par- 
ticularly a parameter called s that is proportional to 

^ '' n 

We here point out that the new discovery [4|] of a neu- 
tron star with a mass equal to 1.97(4) solar masses (a 
convenient unit weighing about 2 x 10"^° kg), confirrning 
previous claims of neutron stars in this mass range 5] is 
so close to the maximum mass that such an object can 
have [6|, by nuclear physics considerations, that it con- 
strains the value that the gravitational constant can take 
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FIG. 1. The gravitational constant remains (so far) a con- 
stant. Newton-Cavendish constant normalized by its accepted 
value 6. 6738(8) A'^(m/fc(7)^. Left point: laboratory on Earth. 
Middle: orbital determinations of binary pulsars. Right: ex- 
istence of a neutron star with mass 1.97(4) solar masses. At 
the intense gravitational field in such neutron star, the grav- 
itational constant cannot exceed 8% of its value on Earth at 
95% confidence level. 



in its interior. The reason is that there is an equilib- 
rium between gravitational attraction and interneutron 
repulsion at short distances, that cannot be maintained 
somewhat above two solar masses, and heavier objects 
collapse into black holes. Thus an increase of the gravi- 
tational constant Gm that diminishes the maximum mass 
attainable, is excluded by the discovery of these super- 
heavy neutron stars (most neutron stars known to date 
had masses near 1.4 times that of the Sun). 

The situation is described in figure [Ij From left to 
right, we show the laboratory value (precision 1 x 10^"*); 
the astronomical value as inferred from double pulsars 
J0737-3039 and Hulse-Taylor PSR B 1913+16, at accel- 
eration near 300 m/s^ (respective precisions 5 x 10~^ 
and 2 x 10~^), and our constraint from the existence of 



a neutron star with mass equal to 1.97(4) solar masses 
(Gjv cannot exceed its earthly value by more than 8% at 
95% confidence level). 

To obtain the bound we employ the equation of hydro- 
static equilibrium of Tolman-Oppenheimer-Volkoff [7J, |8|, 
consequence of Eq. ([2]). This equation governs the vari- 
ation of the pressure P inside a spherically symmetric, 
static star at radial distance r from its center, inside 
which a mass M{r), due to the mass-energy density e(r), 
has accumulated 
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dr 
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This equation can be integrated numerically from the 
inside of the star (r = 0) to the outside by a standard 
Runge-Kutta computer algorithm. The initial condition 
is supplied as a value of the pressure in the star's center, 
and the equation is considered solved at the distance R 
where the pressure drops to zero. R is then interpreted as 
the star radius. The mass function M{r) and total star's 
mass M{R) are obtained by adopting an equation of state 
that relates the total energy density e(r) to the pressure 
P{r). Here is where steady progress in nuclear theory 
allows us to have reasonable confidence in the equation 
of state input, shown as the solid red- line in figure [21 
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FIG. 2. The equation of state (pressure as function of density) 
for pure neutron matter (solid red line), compared with the 
free neutron gas (dashed black line) is much stiffer because of 
repulsive interactions (nuclear matter is barely compressible) 
and an independent low-density effective theory equation of 
state [3]. Also shown is the causality limit, given by the condi- 
tion that sound propagates slower than light, c^ — dP/dp < c. 



The plot also shows, for comparison the equation of 
state for a pure neutron Fermi gas [l3|- This is much 
less "stiff" (lower pressure at given energy density) since 
it does not include the interneutron interaction, that is 
repulsive. Further, we show an independent computation 



of the equation of state at low energy densities that is in 
good agreement @ with our own computation. A slight 
discrepancy can be adscribed to our employing pure neu- 
tron matter for simplicity while those authors are includ- 
ing a small amount of protons in dynamical /3-equilibrium 
with the neutrons. Our improvements concern the use of 
a recently developed chiral effective field theory (EFT) 
for nuclear matter |lll4l4l | based on a new power count- 
ing. The latter is able to single out the set of Feynman 
diagrams, including infinite strings of them, that are re- 
quired to calculate order by order in the chiral expansion 
for (asymmetric) nuclear matter. This is a novelty since 
previous calculations employ standard many-body meth- 
ods not based in the new paradigm of EFT. The chiral 
power counting [ll| considers multi-nucleon forces both 



from pion exchanges as well as from short ranges contri- 
butions. It also takes into account the important infrared 
enhancement affecting nucleon propagators in the multi- 
nucleon reducible loops. This allows to control the size 
of many-body (three-body, four-body, etc.) forces and 
decide, systematically and to the precision desired, what 
medium effects are to be kept at each stage of the calcu- 
lation. 

Further details on the equation of state that we employ 
are documented in Figure [31 that presents the pressure 
and energy densities as functions of the Fermi momen- 
tum for the neutron gas, as well as the speed of sound 
\J dPjdp. It should be noted that, while obtained with 
sophisticated modern effective field theory treatments, 
the equation of state is in broad agreement with vin- 
tage nuclear theory treatments based on phenomenolog- 
ical potentials [15| . 

Once the equation of state has been fixed, and the inte- 
gration of the Tolman-Oppenheimer-Volkoff equilibrium 
equation ([3]) has proceeded, one obtains the standard 
mass-radius plot in figure [4] As is known since the early 
work of Oppenheimer and Volkoff, a pure neutron gas 
supports no star with mass above 0.6-0.7 solar masses, 
providing a check of our computer programme (dashed 
red line). The full calculation including interactions can 
elevate the maximum mass above 2 solar masses. We 
discontinue the solid black line at a point where the ef- 
fective theory breaks down as manifested by reaching the 
causality limit |16| Cg = c (the Fermi momentum at that 
point, about 600 MeV, is also quite high). Stars above 
2.2-2.3 solar masses are not supported. For a star with 
mass about two solar masses, near the maximum possi- 
ble, we give in figure [SI the profiles of pressure and inten- 
sity of gravity from the center of the star. The accelera- 
tion of gravity does not grow uniformly from the center 
to the edge as in Newtonian mechanics (as seen easily 
from Gauss's law) due to the pressure contribution in 
the relativistic expression for the potential $ (given in 
geometrodynamic units c = Gn ~ 1) 
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r{r - 2M(r)) 
The order of magnitude of g in the full calculation can 
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FIG. 3. Equation of state for pure neutron matter in effec- 
tive theory (red solid line) versus the free neutron gas (black 
dashed line). Top: Pressure against Fermi momentum. Mid- 
dle pannel: density as function of the Fermi momentum. Bot- 
tom: speed of sound showing the causality limit where the 
effective theory description is expected to break down. 



be understood from simple Newtonian considerations as 
GnM{R)/R^ at the star's surface, and if corrected by 
the relativistic denominator (1 — 2GnM(R)Ic^)~^ ~ 1.4 
one obtains about g ~ 2 x 10^^ m/s^. [l3| 

Thus, any information gained on Gn extends our 
knowledge several orders of magnitude in field intensity. 
This can be of use to constrain modified theories of grav- 
ity, motivated by string theory and by cosmology [18|, 
that suggest that for largely different values of the cur- 
vature R (or acceleration field g), gravity separates from 
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FIG. 4. Mass-radius plot for the neutron star solutions of the 
Tolman-Oppenheimer-Volkoff equation of hydrostatic equilib- 
rium. As is well known, the free neutron gas equation of state 
cannot reach masses beyond about 0.6 solar masses. However 
the interacting equation of state, being more repulsive, sup- 
ports stars slightly above 2 solar masses against gravitational 
collapse, in agreement with the observation of a star with 
1.97(4) solar masses. 



its Einstcnian formulation. 




FIG. 5. For a neutron star with mass near the maximum al- 
lowed by hydrostatic equilibrium, we show the pressure profile 
and the acceleration of gravity (note its non-Newtonian be- 
haviour due to the relativistic pressure term). 

We complete our analysis by returning to Eq. (j3|) and 
varying Gn- Since reducing it simply delays gravitational 
collapse and eventually allows for arbitrarily heavy stars, 
no constraint is put in smaller-than-physical Gm values, 
as refiected in figure [T] However, increasing Gjq rapidly 
reduces the maximum possible mass of the neutron star. 



To control the systematic uncertainty at very high en- 
ergy densities, where other phenomena might arise (ac- 
tivation of the strangeness degree of freedom, transi- 
tion to a different phase of nuclear matter not acces- 
sible from the nucleon effective theory, etc.), and since 
we are interested in imposing an upper bound on Gat, 
we substitute our equation of state by the stiff-most 
allowed by causality [I3, such that Cs — c, yielding 
P = c^{p ~ Pmax) + ^max, abovc a maximum Fermi mo- 
mentum of either kp = 600 or 450 MeV. Shown in 
figure [5] is the mass/radius plot adopting the first value, 
allowing for the Cavendish constant to vary. From this 
calculation we derive the bound on an 8% variation of 
Gn- 

Should one adopt the second value due to putative er- 
rors that we may have not identified in our equation of 
state at higher energy density, the constrain on Gn is 
somewhat relaxed, but remains meaningful, excluding a 
variation of 19% at the 2a level. 

On the other end, we neglect corrections due to the neu- 
tron star skin containing several atomic sheets [2^ as 
well as ordinary nuclear (not neutron) matter, for it is 
known that its total contribution to the star's mass is 
rather small [i,[2l|. 
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FIG. 6. Mass-radius plot as in figure [4] but varying the grav- 
itational Newton-Cavendish constant. As this grows (or, as 
shown, the ratio of the constant in Earth to that constant at 
high field decreases) the star becomes more prone to gravita- 
tional collapse and thus the maximum reachable mass drops 
below the requisite two solar masses. Thus, increases of Gn 
are now constrained. 

In conclusion, we believe that we have made a rele- 
vant contribution in employing the new two-solar mass 
neutron star to constrain the gravitational constant Gn 
given hard information on the nuclear matter equation 
of state, as opposed to attempting to constrain the equa- 
tion of state and possibly exotic forms of nuclear matter 
about which abundant literature exists |2^i2Q] . 
Other authors p] have already pointed out that nuclear 
physics is precise enough to constrain the radius of the 
star given its mass. The discovery of higher- mass neutron 
stars allows us to establish limits on allowed variations 
of gravity itself in a hitherto unexplored regime. 
Not long ago, a two-solar mass neutron star was thought 
unlikely [27J . Although we now know that nuclear physics 
can accommodate it, the margin is narrow and indeed al- 
lows to constrain gravity. 

Future work may include an examination of modified the- 
ories of gravity (for the case of scalar modifications of the 
action in f{R) theories the Tolman-Oppenheimer-Volkoff 
equations are already available in the literature |28j). 
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